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NOMENCLATURE 


c speed of sound 

D jet diameter 

c t - turbulence model coefficients 

1. f 1 ean and fluctuating components of variable f 

g time factor (Eq. 4 1 ) 

k turbulence kinetic energy, Vi u,Uj 

1 turbulence length scale 

M, turbulent Mach number, Eq. 23. 

Pij stress production, Eq. 6, P = P kk /2 

p pressure 

R radius 

R g gas constant 

Rjj space factor, Eq. 41 . 

R, turbulence Reynolds number 

T temperature 

T,j Lighthill stress tensor 

u,v,w fluctuating velocities in x,y,z directions 

u T friction velocity 

x,y,z Cartesian coordinates 

y + nondimensonal distance, yu T /v 

a factor in Sarkar compressbility correction, Eq. 23 & 43 

P factor in compressible dissiapation model, Eq. 43 

5,j Kronecker delta 

e dissipation rate of turbulence energy 

y ratio of specific heats 

k factor in Zeman compressibility correction, Eq. 43 

|i dynamic viscosity 
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v kinematic viscosity 

v t turbulent viscosity 

0 polar coordinate 

p density 

a turbulent Prandtl number for diffusion of k and e 

x 0 characteristic time daiay 

x T shear stress 

Q source frequency 




SURVEY OF TURBULENCE MODELS FOR THE COMPUTATION OF 
TURBULENT JET FLOW AND NOISE 


Abstract 

The report presents an overview of jet noise computation utilizing the computational fluid 
dynamic solution of the turbulent jet flow field. The jet flow solution obtained with an 
appropriate turbulence model provides the turbulence characteristics needed for the 
computation of jet mixing noise. A brief account of turbulence models that are relevant 
for the jet noise computation is presented. The jet flow solutions that have been directly 
used to calculate jet noise are first reviewed. Then, the turbulent jet flow studies that 
compute the turbulence characteristics that may be used for noise calculations are 
summarized. In particular, flow solutions obtained with the k-e model, algebraic 
Reynolds stress model, and Reynolds stress transport equation model are reviewed. 
Since, the small scale jet mixing noise predictions can be improved by utilizing 
anisotropic turbulence characteristics, turbulence models that can provide the Reynolds 
stress components must now be considered for jet flow computations. In this regard, 
algebraic stress models and Reynolds stress transport models are good candidates. 
Reynolds stress transport models involve more modeling and computational effort and 
time compared to algebraic stress models. Hence, it is recommended that an algebraic 
Reynolds stress model (ASM) be implemented in flow solvers to compute the Reynolds 
stress components. 


I. INTRODUCTION 


The study of the origin of jet noise began around 1950 in response to the then emerging 
need to control the noise of jet propelled aircraft. Lighthill [1,2] proposed the theory of 
aerodynamic sound to describe the mechanism of noise generation from the mixing zone 
of turbulent jets. It is an exact formulation based on the fundamental equations of fluid 
motion. Lighthill’s equation for density fluctuations in a flow is written as, 


dr dX'Xj 

where 

Tjj = pUjUj - Xjj + (p pc‘)5ij 


( 1 ) 


( 2 ) 


is the Lighthill stress tensor. 
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The theory replaces the actual flow by a flow at rest with an acoustic field in which 
waves propagate at constant speed c. The source field for the waves is a quadrupole 
distribution and the strength of the quadrupole in unit volume is given by Lihgthill stress 
tensor, Ty. The double divergence on Tjj indicates that the source is a quadrupole. Thus 
in Lighthill’s analogy the sources move instead of the fluid. Hence if Ty is known 
throughout the real flow field the wave equation (1) can be solved, to evaluate the small 
scale jet mixing noise. 

In the absence of a detailed flow field solution, simple scaling laws derived for the 
turbulent flow were used to estimate the sound radiation from turbulent jets. Such 
estimates showed poor agreement with the data. The variation of sound spectra with 
angle to the jet axis was poorly estimated at moderate and high frequencies. This was 
traced to the neglect of mean flow effects on the radiated field. Lilley [3] formulated the 
jet noise problem in terms of jet noise generation and sound-flow interaction, accounting 
for the effect of refraction and convection. This formulation is used in a majority of 
recent investigations of jet noise. 

The problem of estimating the distribution of Ty throughout the flow field has been the 
subject of numerous investigations. A fully time dependent numerical simulation (Direct 
Numerical Simulation, DNS) of the turbulent jet flow can be used to provide the 
distribution of noise source strength Ty. But such full simulations are still restricted to 
simple flows and low Mach numbers. Colonious et al [4] computed the acoustic field due 
to plane mixing layer, using direct simulation. Before we look for other methods of 
computing Ty, let us look at the terms in Lighthill stress tensor. The first term pu ; Uj is the 
momentum flux per unit volume. The second term -Tij is the viscous stress, which can be 
neglected for high Reynolds number flows. The third term (p - pc 2 )8jj is normally 
considered to be small order compared to pujuj in isentropic flows, where the temperature 
difference between the flow and the ambient is small. So for majority of the flows of 
interest Ty =pujUj. pujUj is the unsteady Reynolds stress. However, the full space-time 
history of Ty can not easily be evaluated for flows of practical interest.. 

The Reynolds stress distribution can be obtained from the solution of Reynolds averaged 
Navier-Stokes (RANS) equations. Substitution of apparent mean (Reynolds) stresses for 
the actual transfer of momentum by the velocity fluctuations increase the number of 
unknowns above the number of equations. The problem then is to supply the information 
missing from the time-averaged equations by formulating a model to describe some or 
all of the six independent Reynolds stresses, -pUjUj. The exact Reynolds stress transport 
equations can be derived from the time dependent Navier-Stokes equations [5]. These 
equations express the conservation of each Reynolds stress as the Navier-Stokes 
equations express the conservation of each component of momentum. In turbulence 
modeling one uses a finite number of Reynolds stress transport equations and supplies 
missing information from experimental (or analytical) results. The time-averaged scalar 
transport equation contains the turbulent heat or mass flux, -puj<|), where (j) is the 
fluctuating scalar quantity. When only time averaged information is available, modeling 
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of the turbulent velocity frequency-wave number spectrum is required to obtain noise 
spectra as a function of directivity angle. 

From a time averaged solution with appropriate turbulence modeling, turbulence length 
and time scales needed for the acoustic solution can be extracted. This approach has been 
adopted by Khavaran et al [6], Bailey et al [7,8], and Khavaran and Krejsa [9] recently to 
compute the sound radiated from turbulent jets. These papers use k-e turbulence models 
and they express the turbulence length and time scales in terms of turbulent kinetic 
energy, k, and its dissipation rate, e. In this report, we will look at the turbulence models 
that provide k and e, for use in noise calculations. 

Several reviews of turbulence models have appeared concentrating on different aspects of 
turbulence modeling [for example, 10 -15]. A recent review by Hanjalic [16] 
summarizes the applications of single point closure methods and discusses possible 
directions for turbulence model improvements. Spieziale [17] discusses mathematical 
aspects of Reynolds stress closure methods. The book (revised, 2 nd Edition) by Wilcox 
[18] contains complete details of turbulence models that are employed in computational 
fluid dynamics computations. 

One recent NASA conference publication [19] presents various turbulence models and 
their applications to subsonic/supersonic flows, wall bounded and free shear flows of 
interest in propulsion. Turbulence models used by various industries and research 
organizations and the results obtained with these models are presented. In another 
NASA/Industry report [20] nozzle flow computational results obtained from five 
different codes (from GE, UTRC, MDC, Boeing, and Glenn Research Center (GRC)) 
with different models were evaluated. The codes were found to produce similar results 
when they used common grids, boundary conditions, and turbulence models. The results 
showed little sensitivity to upstream turbulence levels, but showed strong dependence on 
the choice of turbulence model and the near wall treatment. 

In the present survey, we examine the turbulence models that are relevant for the 
computation of jet flow for the purpose of evaluating sound radiated from turbulent jets. 
First a description of turbulence models which are relevant for computing the noise 
radiated by jets is given. Then the application of the models and their performance in jet 
flows of interest are described. 


II. TURBULENCE MODELS 


The transport equations for the Reynolds stress tensor can be derived from Navier-Stokes 
equations [5]. Since such transport equations contain higher order correlation terms, 
models need to be developed to express them in terms of known or calculable variables. 
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2.1 Reynolds Stress Transport Equation Model 

Turbulence models employing transport equations for ujuj are called second order closure 
smodels. Several closure schemes have been proposed for these equations. The well- 
tested one is that of Launder et al [21]. This model was applied to axisymmetric free 
shear flows by Launder and Morse [22]. The free-shear flow version of the transport 
equation for Reynolds stresses transport equations may be expressed as 


Du t u ( 
Dr 


= P„+<l>, i -£ iJ +D, J 


(3) 


Convection = Production + Pressure strain + Dissipation + Diffusion 

The four terms on the right hand side represent the stress production, pressure-strain 
correlation, viscous dissipation and diffusive transport of uiuj, respectively. 

The pressure-strain correlation is approximated as: 


*u =- c . 


(- 


U i“j ~2 5 u k 



f du. 


du , 


(4) 


dx , dx. 



(5) 


where 



/ 


P o=~ 


du , 




dx. 


du i 

dx. 


( 6 ) 




■a«, 


u i u *^f- +u / u t 


du k 


dx. 


dx. 


1 J 


(7) 
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The coefficients a, (3, and y are related to a quantity C 2 by 
a = (8+c 2 )/l 1; |3 = (8c 2 -2)/1 l;y= (30c 2 -2)/55. 

P= P kk /2; k = u k u k /2. 

In [21], two models were adopted for the diffusive transport of stress, Dq. The simpler 
one proposed by Daly and Harlow [23] was used for axisymmetric thin shear layers by 
Launder and Morse and it is: 


D ” =C 'S- 


dit'U t 
dx t 


( 8 ) 


Closure of Launder et al model [21] is completed through the following equation for the 
turbulence dissipation rate, e, of turbulence energy. 


Dz 


Pz 


f, 


~ c ’el , c e2 , + c e-> 
Dt k k ax/. 


-dz 

-UfUk — 

z ox/ 


(9) 


The model contains six coefficients and their values are [22]: 


Ci c 2 c s C E | c e2 c e 
1.5 0.4 0.22 1.45 1.9 0.15 


Convective transport and production terms are exact whereas the diffusion, pressure- 
strain, and viscous dissipation terms have been modeled. The diffusion fluxes of UjUj 
have been expressed by simple gradient diffusion models. The most important 
assumption concerns pressure-strain terms, since for shear stresses these are the main 
terms to balance the production of these quantities. The pressure strain model consists of 
two parts. The first one represents the interaction of fluctuating components only, and 
the second, the interaction of mean strain and fluctuating quantities: <j),j = 1 + 


f. 




\ 


«,■« j -3 s u k 


( 10 ) 


(/> ij — —CL 


P --S„P 

* 3 1 


-P 


( 






^ du. 


<)u , 


(ID 


dx . dx, 


Several versions of pressure-strain model have been proposed to correctly predict the 
experimentally observed results. To account for the wall damping effects a wall 
correction must be introduced in the pressure-strain model. Launder et al [21] make the 
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empirical constants in the pressur-strain model a function of the relative distance from the 
wall, 1/y k ~/(ey). Because of the complexity and the large amount of computational 
effort involved, the model has not been widely used as one would like it to be. 

2.2 Algebraic Stress Model 

In Reynolds stress models, there are differential equations for each component of TijUj in 
addition to an e equation. To reduce computational effort algebraic relations have been 
proposed by Rodi [24] for calculating the Reynolds stresses. This done by assuming that 
the net transport of u,uj is proportional to the net transport of k multiplied by the factor 

UjU/k. 

Rodi uses a simpler model for pressure-strain relation than that presented in Eq. (4) and it 
is given by 


< p = —c. — 

M/ 1 


u i u J 


-« 


! 2 ^ 
P -- 8 P 

{ " 3 " 


( 12 ) 


with a = 0.4 and he writes the transport equation for turbulent energy, k as 


Dk d 

= c. 


( , 


Dr ' ?)x k 


- dk 
-u k u,— 
£ ox, 


-u k u, 


3^ 

dx k 


Dk 


P=P»/2 


(13) 


As mentioned above, to obtain an algebraic expression for UjUj, the following 
approximation is employed: 


Dum 


u.u , 


D =- 4 - 2 . 

Dt 11 k 


Dk 

Dt 


- D, 


u,u , 


( P ~£ ) 


(14) 


where D,j is defined in Eq. (8), Dk and P in Eq. (13). Incorporation of Eq. (14) into the 
UjUj equation (3) yields the desired algebraic expression for UjUj : 


u,u , = k 


0 . P„/£ — S u P/e 

3 ,i+ T 1+ up N 

^ I 


(15) 
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Now we have a set of algebraic expressions for the stresses "uiuj, in terms of the mean 
strain rate, turbulent kinetic energy k, and its dissipation rate e, and the stresses 
themselves. As in Launder et al model [21], closure is completed by an equation for the 
dissipation rate of turbulence energy, e. 

The algebraic stress model provides a mechanism by which anisotropic turbulence 
distribution can be computed without the large amount of computational effort required 
for the Reynolds stress transport equation model discussed above. All the effects that 
enter the transport equations for ujuj through the source terms for example, body force 
effects (buoyancy, rotation, and streamline curvature), non-isotropic strain field and wall 
damping influence can be incorporated into algebraic stress models. Algebraic stress 
models therefore also simulate many of the flow phenomena that were described 
successfully by Reynolds stress transport equation models. 

2.3 k-e Model 

The k-e model is the most often used model in present day engineering computations. 
The model was developed by Launder and Spalding [25,26] and Hanjalic and Launder 
[27]. In this model closure is achieved by relating the Reynolds stress to the mean strain 
rate through the Boussinesq approximation 


- pUjUj = p, 


( du. ] 


-+- 


dx , dx, 


-—S'jpk 
3 1 


(16) 


The effective turbulent viscosity, p, is defined in terms of a characteristic length and 
velocity. If the length scale is taken as the turbulent length scale, k V2 /e, and the velocity 
scale is approximated as Vk, then p, can be expressed as 

p,=Cnpk 2 /e (17) 

c M is a constant. The individual ujuj is related to the single velocity scale V k. For 
isotropic turbulence ujiij =2/3 5jjk. In k-e model one solves two separate modeled 
transport equations, one for turbulent kinetic energy and the other for its dissipation rate. 
The modeled equations for k and e as described in Reference 26 are: 

(a) Kinetic energy equation 


Dk i a 

p f dk 

+ Hl 

f du, du k ^ 
1 

diif 

Cl 

1 

1 +-* 
|q 

_o* dx k J 

P 

dxi dx: 

k 1 ) 

dx k 


( 18 ) 


(b) Kinetic energy dissipation rate equation 
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(19) 


De 

~Dt 


1 8 

P, de 


'du, du, ^ 

1 -A- K 

8m, 

P dx k 

1 

P k 

1 

dx, dx, , 

A ' J 

dx k 


C c2 


£' 


k 


The constants assume the approximate values of = 0.09, c e i = 1 .44, c e 2 = 1 .92, ctk = 1 .0, 
and o f = 1.3. These constants were obtained by comparison of model predictions with 
the experimental data on equilibrium boundary layers and decay of isotropic turbulence. 


III. MODIFICATIONS TO k-e MODEL 


The standard k-e model has been modified to account for observed discrepancies between 
the model prediction and the experimental results. Here we consider first two such 
modifications relevant for the computation of jet flows to account for the spreading rate 
of circular jets and the spreading rate of high-speed jets. Then we discuss an anisotropic 
k-e model, low Reynolds number and near-wall models, and multiple-scale models. 

3.1 Vortex stretching dependent dissipation rate 

It was found early on that while the standard k-e model predicts the plane jet flow 
correctly, it overestimates the spreading rate of circular jets. Pope [28] suggested that the 
stretching of vortex tubes by the mean flow has significant influence on the process of 
turbulence scale reduction. In axisymmetric jets, as the jet spreads rings of vorticity are 
stretched. This causes the effective viscosity and hence the spreading rate to be lower in 
the circular jet. Pope incorporated this aspect in the standard k-e model by modifying the 
dissipation rate, e, equation. The modified form of the dissipation equation proposed by 
him is: 


De 

Dt 


1 8 

P, de 

, c eI /i, £ 

8m 8m, ^ 

1 . A 

8m, 

P dx k 

a e dx k _ 

P k 

1 

dx, dx, 

V * 1 ) 

dx k 



Cc)X 


( 20 ) 


Where % = cOjjttySjj 


1 * 

S ‘ J ~ 2 e 


fa«, + a«, A 


v 


dXj dx t j 


1 k 

0)0 ~ 2 e 


f du t du j A 


dxj dx, j 


( 21 ) 


( 22 ) 
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and Cg 3 = 0.79. 


3.2 Compressibility Correction 

The standard k-e model when used to predict the development of high-speed shear layers 
and jets, it was found that the growth rate did not compare well with the measurements. 
In these flows, the experiments showed that the growth rate of high-speed shear layers 
reduces with increase in convective Mach number [29]. The growth rate of shear layers 
is dependent on the growth rate of instability waves at these speeds. At high speeds, the 
reduction in the instability wave growth rate reduces turbulent mixing. Sarkar et al [30], 
Sarkar and Lashhmanan [31], and Sarkar [32], developed an addtitonal factor to be added 
to the standard k-e model to account for the compressibility effects. The form of the 
additional factor was found by an asymptotic analysis of the compressible Navier-Stokes 
equations. The suggested modification is 

e= ed 1+ocM, 2 ) (23) 

where M, 2 = 2k/(yR gT ) and a is a constant set equal to 1 .0 and R g is the gas constant. 

The factor e= e s (l+aM, 2 ) corresponds to the contribution due to the incompressible and 
compressible dissipation rates, e* referring to the standard value and M t is the turbulent 
Mach number. This term is added to the turbulent kinetic energy equation of the standard 
k-e model. The equation now reads as 


Dk _ 1 3 

p, dk 


( dui t 
1 

g| 

1 

“O 

a* 

_O k dx k _ 

P 

dx k dx: 

K 1 ) 


-efl + aM, 2 ) 


(24) 


3.3 Anisotropic k-e Model 

The standard k-e model assumes an isotropic eddy viscosity relationship for the Reynods 
stress tensor. Reynolds stress models discussed above can predict the observed 
anisotropy in normal stresses. Anisotropic k-e models based on anisotropic eddy 
diffusivitites have been proposed [33-38]. The anisotropic model proposed by Myong 
and Kasagi [35], is valid up to the wall. In this model, the deviations from isotropic 
Reynolds stresses are given by a function of nonlinear quadratic terms of mean velocity 
gradients and that of anisotropic diffusion terms of turbulent kinetic energy. The normal 
Reynolds stresses are algebraically calculated. The expression for Reynolds stress is 
given as: 


Dk 

Dt 


3x 


3 dk 

(V+V, / <T t ) u.u 


dx , 


du 

dx. 


-£ 


(25) 
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f 2 = (l-^exp[-(-§-) 2 ])[l-exp(-^— )] 2 ( 31 ) 

9 6 5 

W = —1.5 — 0.75(5 S + 8 8 ) + 2 (8 8 +8 8 ) + 8 8 8 +8 8 8 (32) 

R,=k 2 /ve ; 8k=1.4, 5,= 1.3, c e! =1.4, c e2 =1.8, and c M =0.09. 

(The indices n and m denote the wall normal and streamwise coordinates respectively). 

The mean velocity, the turbulent kinetic energy and its dissipation rate are not influenced 
by the normal stress anisotropy. The transport equations to be solved are similar to those 
of isotropic k-e model. 

Myong and Kasagi [35] showed that their anisotropic model predicts correctly the 
dependence of each normal component of Reynods stress correctly, u a y, v a y 2 , and w 
a y [Fig. 1]. For the flow over a flat plate, the model predicts the wall-limiting behavior 
that is in good agreement with the data [Fig. 2], The predicted Reynolds stress 
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Fig. 1 Wall-limiting behavior of three normal Reynolds stresses 
- after Myong and Kasagi [35] 



Fig.2 Distributions of normal Reynolds stresses in a turbulent 
boundary layer - after Myong and Kasagi [35] 
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components in the entire region were also found to agree fairly well with the 
experimental results. 

3.4 Low Reynolds Number and Near Wall k-e Model 


Jones and Launder [36], extended the k-e model to model low Reynolds number flows so 
that the turbulence model equation can be valid throughout the laminar, transition, and 
fully turbulent regions. In this version of the model k and e are determined from the 
following equations: 


Dk 

~Dt 


f Vt ) 

dk 

+ Hl 

dui ( du k 1 

Bu: 

fa* 1 ' 2 ) 

LI ** J 

dx k _ 

p 

, dx k dx, . 

LV 

dx k 

Bx : 

1 J 


(33) 


De 

Dt 


’fu+iL) 

de 

, c e\ Pt e 

( duj | Bu k ^ 

d“i , e 2 , Q v Pr 

( ^2- ^ 
d Uj 


dx k 

P k 

{ dx k dx, } 

dx k k p 

dxjdx: 

) 


(34) 


(X, is the turbulent viscosity defined, for the standard k-e model, in Eq. (17). In this 
model, c M and c E 2 vary with turbulence Reynolds number, R t . 


R, = pk 2 /pe 


(35) 


c M = c MS exp[-2.5/( l+R,/50)] (36) 

Cg 2 = c e2s [1-0.3 exp(-R, 2 )] (37) 


Subscript s refers to the standard model values. We note here, that the laminar diffusive 
transport becomes of increasing importance as the wall is approached and the extra 
destruction terms included are of some significance in the viscous and transitional 
regions. The term. 


2.0 




dx,dx. 


in the £ equation produces satisfactory variation of k with distance from the wall. In the 
computations £ is set to zero at the wall and an extra term. 


f 

-2v 

V 


3A 1/2 

d Xj 


Y 

y 


is introduced to the k equation. This extra term is exactly equal to the energy dissipation 
rate in the neighborhood of the wall. 
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Turbulence models for near wall and low Reynolds number flows were reviewed by Patel 
et al [37]. Eight different models ( all based on k-e model except one) were considered 
and their performance in predicting turbulent boundary layers with and without pressure 
gradient (favorable/adverse) was examined. The model of Launder and Sharma [38] and 
that of Chien [39]. both based on Jones and Launder model described above appear to 
perform well in majority of the test cases studied by Patel et al [37]. 

The model of Chien [39] is claimed to perform better than that of Jones and Launder is 
briefly described here. Though the model is based on Jones and Launder model, the 
presence of solid wall is handled differently. An additional term, representing the finite 
dissipation rate at the wall, is added to balance the molecular diffusion term. The 
dissipation term in the kinetic energy equation is given by e + (2vk/y“) for finite values of 
y, distance from the wall. The turbulent kinetic energy equation takes the form 


Dk 

Dt 


d_ 

dy 


(v + v f ) 


dk 

dy 


du 2 

+v 'V 


2 vk 
-e - 

V 


(38) 


The term 


2vk 

i 

y 

is the term added to produce correct behavior of turbulent energy k in the near wall 
region, v is the kinematic viscosity. The turbulent viscosity v t is modified to reflect the 
wall damping effect. 


v, =c — (l-exp(-c 3 w*;y/v) 


(39) 


c 3 is a constant, u* is the friction velocity. The turbulent dissipation rate equation 
suggested by Chien reads as 


De _ 9 

Dt ~ dy 


, v de 
(v + — )— ■ 

c„ dy 


£ ,9m. -> £ 

+e -rV* 


c El f£ + 


2 vk exp(-c 4 w y / v ) 


(40) 


where f = 1-0.222 exp[-(R,/6) 2 ], c 4 is a constant. c 3 = 0.01 15 and c 4 = 0.5 were used by 
Chien. 
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3.5 Multiple-Scale Model 


The turbulence models discussed above are based on the assumption that in all flow 
situations turbulence has a spectrum of universal form which can be characterized by the 
scale of the energy containing range. Difficulties arise when the spectrum is not an 
equilibrium one or when the flow exhibits distinctly different ranges of scales. A two- 
scales model was proposed by Hanjalic et al [40] . They split the spectrum into a large 
scale part and a small scale part with different time scales for energy transfer into the 
large scale part and transfer from large scale to small scale part. 

The turbulence spectrum consists of independent production, inertial, and dissipation 
ranges. Ki denotes the wave number above which a significant mean strain production 
occurs while Ki is the largest wave number at which viscous dissipation of turbulence is 
unimportant (Fig. 3). Energy leaves the first region (production) at a rate e p and enters 
the high wave number or dissipation region at a rate £,. Between the two regions, 
occupying the intermediate range of wave numbers is the transfer region, across which a 
representative spectral energy transfer rate e T is assumed. This simplified energy 
spectrum is the basis of the model of Hanjalic et al. The total turbulence energy k is 
assumed to be divided between production range k p and the transfer range k T At high 
Reynolds numbers there is negligible kinetic energy in the dissipation range. The 
transport equations for k p , k T . £ p , and e T are formulated. Thus there are two k and two £ 
equations in this model and two sets of constants which are determined from 
experiments. 



Fig. 3 The spectral division for multiple scale model - after Hanjalic et al [40] 
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Modified versions of the above two-scale model have been formulated by Kim & Chen 
[41] and Chen [42], In the model used by Duncan et al [43], the model coefficients were 
made dynamically dependent on the partitioning of the energy spectrum. Ko and Rhode 
[44] developed a new multi-scale k-e turbulence model, which incorporated a new way of 
evaluating source/sink coefficient functions. Though these models are attractive from a 
theoretical viewpoint, their use to flows of engineering interest is hampered by the 
number of constants needed to be calibrated with these models. 

Next, the application of the turbulence models to the prediction of jet nose shall be 
discussed. 


IV. TURBULENCE MODELS AS APPLIED TO JET NOISE PREDICTION 


4.1 k-e Model Computations for Jet Noise Prediction 

The quadrupole source term (unsteady Reynolds stress) that appears in Lighthill’s 
equation has to be evaluated to compute the jet noise. In the absence of detailed time 
dependent flow information, one uses the mean flow information from a simplified 
turbulent flow model such as that of Reichardf s[45]. Suggestions were made that with 
the advances in computational fluid dynamics (CFD), the source terms can be computed 
more accurately from the solution of Reynolds averaged Navier-Stokes equations using a 
k-e model[46], Khavaran et al [6] were the first to carry out such a source computation 
and use the source characteristics for the computation of jet noise. They considered a 
convergent-divergent nozzle geometry. The flow solution was obtained using an 
axisymmetric version of PARC code [47] with Chiens’s k-e model [39]. They showed 
good agreements of the CFD results with the data. The computed turbulence intensity 
contours in the flow field are shown in Fig. 4. Comparisons of the computed turbulence 
intensities with the data and the Reichardt’s solution are shown in Fig. 5. 



Fig. 4 Tubulent intensity contours in a round jet — after Khavaran et al [6] 
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Fig. 5 Comparison of predicted and experimental turbulent intensities 
(a) Lip-line and (b) X/D - 8.2 1 - after Khavam et al [6] 


The time averaged flow information is used to compute the sound field: The noise source 
strength in a turbulent flow is characterized by a two-point time delayed fourth order 
velocity correlation tensor. The fourth order correlation is expressed as a linear 
combination of second order correlations. Then, the two point velocity correlations are 
written in terms of separable space/time factors as suggested by Ribner [48], 

UjUj = Rjj(% )g( T) (41) 

The space factor Rjj(^) is expressed as a function of turbulence intensity and the 
longitudinal macroscale of turbulence. The time correlation is expressed in terms of the 
characteristic time delay, t 0 , which is proportional to the inverse of mean shear and is 
related to turbulence kinetic energy k (k = uiuj/2) and its dissipation rate e as t 0 = k/e. 
Thus, the noise source strengths can be expressed in terms of length and time scales 
extracted from time averaged solutions. The corresponding spectrum can then be 
evaluated by a using a Fourier transform on the time delay of correlation. 

Khavaran et al evaluated the contribution to self noise for various source strength 
components using Ribner’ s formulation [48], The contribution to the acoustic pressure, 
p(R,9,£2), due to each quadrupole source may be expressed as 


P 2 (R,0,Q) oc 



(42) 


where R is the radius, 0 is the angle with respect to jet axis and Q the source frequency. 
It is seen that for accurate prediction of acoustic pressure, the turbulent kinetic energy 
and its dissipation rate £ need to be computed accurately by the flow solution (x 0 = k/e). 
Further details of the noise computation can be found in [6]. Comparisons of the overall 
sound pressure level directivity show good agreement with the data and the results of 
Reichardt method (Fig. 6). The spectral components of noise, based on the one-third 
octave band, are shown in Fig. 7. 


16 






100 1 1 1 1 1 1 1 * 

0 20 40 60 60 100 120 140 160 160 

Anglo from axis, 0 


Fig. 6 Comparison of overall sound pressure level directivity with data 
on a 40-ft radius - after Khavaran et al [6]. 
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Fig. 7 Comparison of the spectral components of noise with data based on one~third octave 
center frequency (band number 1 1 corres >onds to 50 Hz). - after khavaran et al [6] 
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Bailly et al [7,8] used a similar approach and computed the jet noise using the source 
strength obtained from a CFD solution using k-e model. They computed far-field levels 
using Ribner’s model and also that of Goldstein and Rosenbaum [49] who modified the 
Ribner model by introducing an anisotropic description of turbulent field. They found 
that Goldstein and Rosenbaum model produced better agreement with the far-field sound 
pressure data. 

For supersonic jets, Bailly et al used Ffowcs-Williams and Maidanic [50] formulation of 
Lighthill source term to account for Mach wave radiation which is one of the main noise 
source in supersonic jets. They applied this model when the local convective Mach 
number is supersonic. Using a combination of Goldstein and Rosenbaum model and 
Ffowics-Williams and Maidanic model, they were able to compute the far-field acoustic 
pressure as a function of jet Mach number, from low subsonic to high supersonic jets. 
The computed far-field levels were found to be in good agreement with the data (Fig. 8). 
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Fig. 8 Total acoustic power as a function of jet Mach number: comparison 
Of hybrid noise model to experimental data - after Bailly et al [7J 
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However, the description of the two-point velocity correlation needs to be further 
examined to improve the source definition and better establish the proportionality factors 
that arise from the use of time averaged quantities. 

Measured turbulence intensities in nozzle flows show considerable anisotropy among 
turbulence intensity components. Longitudinal component (U|) is significantly different 
from the transverse components (u 2 and u 2 ). The transverse components (u 2 and u 2 ) are 
nearly equal. To account for the observed anisotropy Khavaran and Krejsa [9] proposed 
the use of an axisymmetric turbulence model for jet noise computation. However, since 
the CFD solutions were obtained with a k-£ model (isotropic turbulence assumed) they 
investigate the influence of anisotropy by varying the ratio of transverse to longitudinal 
length scales, Ai, and the ratio of intensities, A u = u 2 2 /u| 2 . They demonstrated the effects 
of the parameters on the noise directivity of a splitter nozzle flow. The predicted noise 
directivity shows good agreement with the data when the parameter A| = 0.5 and A u = 0.6 
(Fig. 9). The effects of anisotropy parameters are summarized in Fig. 10. An increase in 
anisotropy tends to increase the sound pressure level. The use of axisymmetric 
turbulence model for noise computation would improve the noise predictive capability of 
jets of practical interest. 



Fig. 9 Predicted sound pressure level directivity for the splitter nozle on a SO ft arc. 
Anisotropy parameters are: A “ 0.50 and UjVui 2 “ 0.60. — after Khavam and Krejsa [9]. 


19 





Fig. 10 Sound pressure level directivity vs. non-isotropy parameters A and 
u 2 J /ui . - after khavaran and Krejsa [9] 


Since, distributions of U| 2 and U 2 2 are needed for the prediction of jet noise, ways to 
obtain these components should be explored. The possibilities exist to use the Reynolds 
stress transport equation model or the algebraic stress model (ASM) described in section 
2 . 

4.2 Other k-e Model Predictions 

Here, some other applications of turbulence models to jet flow predictions that produce 
reasonable solutions that may be used for jet noise prediction are discussed. 

Numerous k-e model predictions have been carried out for jet flows. But these 
predictions were mainly intended to study the flow field characteristics and they have not 
been used for the purpose of noise prediction. Some of them are reviewed here, as they 
hold promise for noise predictions. 

One of the most exhaustive applications of k-e model for jet flows encompassing 
subsonic, supersonic, cold and hot jet flows is that of Thies and Tam [51]. The jet Mach 
number varied from 0.4 to 2.2 and the ratio of jet reservoir temperature to ambient 
temperature varied from 1.0 (cold jet) to 4.0. They demonstrated that if the original 
constants of the k-e model are replaced by a new set of constants (established 
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empirically), the calculated jet mean velocity profiles agreed with the data for a wide 
range of jet flows. They included the correction term for vortex stretching due to Pope 
[28] and compressibility correction due to Sarkar [32], but with new empirically 
established constants. Their choices of model constants are: 

Cjj C, i C e 2 Qk &£ CL 

0.0874 1.4 2.02 0.822 0.324 0.377 0.518 

Note that the factor associated with the vortex stretching term, c e 3 . and the factor, a, 
associated with the compressibility correction term are also modified. The parabolized 
equations, in the Favre-averaged form, are solved using an accurate dispersion-relation- 
preserving (DRP) numerical scheme. In all the cases the computation started from the 
nozzle exit, with initial conditions derived analytically or from the data. The predicted 
mean velocities agreed well with the data as shown in Fig. 1 1 for heated jets. 

Dash et al [52-54] in a series of papers have explored different formulations of k-£ model 
and its various combinations for jet flow predictions. A k-e model with modified 
compressible dissipation factor (due to Sarkar [32] and Zeman [55]) and with Pope 
correction factor was found to yield reasonably good predictions over a range of jet flow 
conditions. They expressed the compressible dissipation as 


e t = e s \piM; + jSM, 4 ] (43) 

where a = 1 (same as Sarkar) and M, = M, - X. 

X = 0. 1 (same as Zeman) 
and P= 60, to fit LaRC data best. 

An example of their predictions of centerline velocity (Fig. 12) and temperature (Fig. 13) 
for different jet exit temperatures of Seiner’s [56] jet are shown. The trends are predicted 
reasonably well. 

The use of compressible dissipation factor for supersonic jet flow predictions was also 
studied by Balakrishnan et al [57], They found that with the compressible dissipation 
correction the reduced spreading rate of supersonic jets was successfully predicted (Fig. 
14). The prediction of pressure distribution in an under-expanded jet with and without 
compressibility correction is shown in Fig. 15. The improvements observed due to 
compressibility correction factor in predicting turbulence intensities in an under expanded 
jet are shown in Fig. 16. for two numerical algorithms. 
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Fig. 1 1 Comparisons of computed and measured center line velocity and half-velocity 
point distributions for Seiner ei al’s [53] Mach 2.0, axisymmetric jet o, asured centerline 
velocity and II, measured half-velocity point. - after Thies and Tam [51] 
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Fig. 13 Axis temperature decay of the Seiner Mach 2 jet cases into still air using the jet 
Modified keCD model with the Pope centerline correction. - after Dash et al [52] 
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Fig. 14 Velocity profile at X/D = 13.5 for a supersonic jet - after Baiakrishnan and 
Abdol-Hamid [57]. 
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Fig. 15 Centerline pressure distribution for 
an underexpanded supersonic jet. - after 
Baiakrishnan and Abdol-Hamid (57) 
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Fig. 16 Lipline turbulence intensity for 
an underexpanded supersonic jet. — after 
Baiakrishnan and Abdol-Hamid [57]. 
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Renormalization group (RNG) theory [58] based k-e model was incorporated into 
NPARC code by Papp and Ghia [59], They compared the solutions obtained with the 
RNG based k-e model with that of Chien's model [39]. They found that while RNG 
model produced slightly better results, the Chien’s model exhibited numerical stability 
problems for certain grid resolutions. RNG method appeared to be more robust. 

Tubofan engine exhaust nozzle flows [60] and lobed mixer nozzle flows [61] have been 
successfully calculated with standard k-e model. 

4.3 Algebraic Stress Model (ASM) Prediction 

Now with the ability to use axisymmetric turbulence characteristics to compute noise 
established [9] it is imperative that the Reynolds stress components should be calculated 
accurately. Algebraic stress model provides that possibility. An example of ASM model 
prediction of a coaxial jet flow [62] is shown in Fig. 17 and 18. The figures show u and 
v profile variations with axial distance respectively. The ratio of turbulence intensities, 
A u = U 2 2 /ui 2 and the ratio of transverse to longitudinal length scales can now be calculated 
accurately for use in the axisymmetric turbulence model for noise calculation [9]. 

4.4 Reynolds Stress Transport Equations Model Prediction 

This model can be used to compute the Reynolds stress components accurately. This 
model requires modeling of higher order correlations and additional computational efforts 
to solve the transport equation for each component of the Reynolds stress. For those 
reasons this model is not as widely used as one would like to be. Figure 19 shows the 
three components of normal stress profiles in a round jet, computed using the Reynolds 
stress transport equation model [22]. 


V. LOCATION OF INLET BOUNDARY AND BOUNDARY CONDITIONS 


The inlet boundary condition specification plays a crucial role in turbulence model 
predictions [63,64]. One can locate the inlet boundary for the jet flow calculations at the 
jet exit and specify the inlet conditions there, as for example was done by Thies and Tam 
[51]. Such a choice, avoids the complexities associated with the boundary region such as 
the near wall modeling and wall limiting behavior of turbulence quantities. However, 
formodern complex and novel designs such as those with internal mixing device, ejector, 
etc. the jet exit flow conditions can not easily be generated accurately (or known from 
measurements). In such complex nozzle flow cases, the nozzle internal flow has to be 
computed with appropriate near-wall modeling in the nozzle-mixer-elector flow path. 
Chien’s model [39] has been used successfully [6,9] for several jet flow configurations. 
But, it is essential to incorporate correct wall limiting behavior turbulence quantities by 
employing a model such as that of Myong and Kasagi [35], The use algebraic Reynolds 
stress model (ASM) will provide an accurate description of Reynolds stress components 
[62], 
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Fig. 17 Predicted and measure profiles of RMS axial velocity, u , for coaxial 
Jets in ambient air. - after Srinivasan et al [62]. 
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Fig. 18 Predicted and measured profiles of fluctuating radial velocity component, v , in 
Coaxial jets. - after Srinivasan et al [62]. 
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Fig. 19 Normal stress profiles in round jets , — usual thin shear flow form, 

including secondary production terms. - after Launder and Morse (22]. 

VI. NUMERICAL SOLUTION ALGORITHM AND TURBULENCE MODEL 


Studies have shown that the same turbulence model incorporated into different codes 
produce different turbulence characteristics [57,59,65,66]. This may arise due to several 
factors such as the numerical solution algorithm, grid dependence, turbulence model 
methodology and implementation, and near-wall model. Flow solvers and turbulence 
models need careful bench mark testing for jet flow computations so that they can be 
used with confidence for acoustic assessment of new nozzle designs. 


VII. CONCLUDING REMARKS 


A brief account of turbulence models that are relevant to provide turbulence 
characteristics needed for jet mixing noise calculations is presented. Length and time 
scales should be predicted accurately to estimate the sound pressure levels correctly. The 
use of compressibility correction due to Sarkar results in correct spreading rates in 
supersonic jets. For axisymmetric configurations, vortex stretching parameter correction 
due to Pope provides the correct jet spreading rate. It is recommended that a near-wall 
model that produces correct wall-limiting behavior of Reynolds stress components be 
used. Anisotropic turbulence information should be incorporated in the small scale 
mixing noise calculation to improve the far-filed noise level estimates and spectral 
distribution. 

Jet flow computations that present the components of Reynolds stress are scarce (as 
indicated by sections 4.3 and 4.4). It is perhaps due to the fact there was no immediate 
use for them. Moreover models such as algebraic stress models and Reynolds stress 
transport models were mostly used for complex flows such as non-circular duct flows, 
curved flows, flows with large separated regions, etc. Recently, it has been shown that a 
knowledge of the magnitudes of the Reynolds stress components is essential for accurate 
evaluation of jet noise levels [7,9]. Turbulence models that can provide the distribution 
of Reynolds stress components must now be considered for jet flow computations. In 
this regard, algebraic stress models and Reynolds stress transport models are good 
candidates. Reynolds stress transport models involve substantially more modeling, and 
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computational effort and time compared to algebraic stress models (section 2.1). Hence, 
it is recommended that an algebraic Reynolds stress model be implemented in the flow 
solvers (such as NPARC code) and validated. Anisotropic turbulence characteristics 
obtained using such a turbulence model would substantially improve the confidence 
levels in jet mixing noise predictions. 
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